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ABSTRACT 

We calculate the advection/diffusion of the large-scale magnetic field threading an advection dominated 
accretion flow (ADAF), and find that the magnetic field can be dragged inward by the accretion flow efficiently, 
if the magnetic Prandtl number V m = rj/v ~ 1. This is due to the large radial velocity of the ADAF. It is found 
that the magnetic pressure can be as high as ~ 50% of the gas pressure in the inner region of the ADAF close 
to the black hole horizon, even if the external imposed homogeneous vertical field strength is < 5% of the gas 
pressure at the outer radius of the ADAF, which is caused by the gas in the ADAF plunging rapidly to the black 
hole within the marginal stable circular orbit. In the inner region of the ADAF, the accretion flow is significantly 
pressured in the vertical direction by the magnetic fields, and therefore its gas pressure can be two orders of 
magnitude higher than that in the ADAF without magnetic fields. This means that the magnetic field strength 
near the black hole is underestimated by assuming equipartition between magnetic and gas pressure with the 
conventional ADAF model. Our results show that the magnetic field strength of the flow near the black hole 
horizon can be more than one order of magnitude higher than that in the ADAF at ~ 3R g {R„ = 2GM/ c 2 ), which 
implies the Blandford-Znajek mechanism could be more important than the Blandford-Payne mechanism for 
ADAFs. We find that the accretion flow is decelerated near the black hole by the magnetic field when the 
external imposed field is strong enough or the gas pressure of the flow is low at the outer radius, or both. This 
corresponds to a critical accretion rate, below which the accretion flow will be arrested by the magnetic field 
near the black hole for a given external imposed field. In this case, the gas may accrete as magnetically confined 
blobs diffusing through field lines in the region very close to the black hole horizon, which is similar to those 
in compact stars. Our calculations are also valid for the case that the inner ADAF connects to the outer cold 
thin disk at a certain radius. In this case, the advection of the external fields is quite inefficient in the outer thin 
disk due to its low radial velocity, and the field lines thread the disk almost vertically, while these field lines 
can be efficiently dragged inward by the radial motion of the inner ADAF. 
Subject headings: accretion, accretion disks, galaxies: jets, magnetic fields, galaxies: active 



1. INTRODUCTION 

The winds driven from the accretion disk through the 
magnetic field lines threading the disk have been consid- 
ered as promising explanations for jets/outflows observed 
in different types of the sources, such as, active galac- 
tic nuclei (AGNs), X-ray binaries , and young stellar ob- 
jects (see reviews in ISoruE [19961 IKonigl & Pudritz 2000; 
iPudritz et all 120071 iSpruitl 120101) . In this model, the or- 
dered magnetic field co-rotates with the gases in the disk, 
and the jets are powered by the gravitation energy released 
by accretion of the gases through the ordered field thread- 
ing the disk (Blandford & Pavnelll982l) . The ordered large- 
scale magnetic field threading the disk is the crucial in- 
gredient in this model. In most of the previous works, 
the strength of the magnetic field is simply assumed to 
scale with the gas/radiation pressure of the accretion disk 
(e.g..lModerski & Sikoralll996l iGhosh & Abramowiczll 19971 
Livio et alJl999tlArmitage & N ataraian 1999; Nemme n et al.l 
20071: IWu&Caol 120081 ; iMcNamara "et all I2009T It 
was suggested that, the external generated large-scale 
poloidal field (e.g., the field originates from the inter- 
stellar medium) would be dragged inward by the ac- 
cretion plasma, while the field will diffuse outward at 
the same time ( Bisnovatyi-Kogan & Ruzmaikin 1974 . 1 1976 



van B allegooiienl 11989c iLubow et alj Il994t lOgilvie & Livic. 
20011) . In this case, the final steady magnetic field configura- 



tion can be derived, in which the inward advection of the field 



lines is balanced by th e outward movemen t of field lines due 
to magnetic diffusion (Lubo w~et alJI 1994b . This means that 
the magnetic field configuration is dominantly determined by 
the radial velocity of the disk and magnetic diffusivity. As the 
radial velocity of the accretion disk is roughly proportional to 
the kinematic viscosity v, the magnetic field configuration is 
sensitive to the magnetic Prandtl number V m —Tj/v, where 77 
is magne ti c dif fusivity. 

Parkei] (119791) argued that v ~ 77 ~ l\\ (I is the largest 
eddy size, and v t is turnover velocity), and Vm = rj/u ~ 
1, is expected in isotropic turbulence. This issue was 
explored by some different authors with numerical simu- 
lations (e.g.. lYousef et alJ 120031: iLesur & Longarettil 12009b 
iFromang & Stond l2009t iGuan & Gammid 120091 which all 
suggest that the magnetic Prandtl number should be around 
unity. A suitable magnetic field configuration is crucial for 
launching a jet from the accretion disk. More specifically, the 
angle of field lines inclined to the midplane of the disk being 
less than ~ 6 0° is required for launching jets from a K eple- 
rian cold disk (lBlandford & Paynell982llCao & Spruit! 19941) . 
This critical angle could be larger than 60° for t he accretio n 
disk surrounding a rapidly spinning black hole (lCaolll997l) . 
which indicates that the spin of black hole may h elp launch- 
ing jets centrifugally by c old magnetized disk s (lCaolll997l 
Sado wski & SikordBOlOl) . ILubow et all (11994 explored the 
final steady magnetic field configuration with the balance be- 
tween advection and diffusion of the large-scale magnetic 
field, and found that significant inward dragging of fields oc- 
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curs only if V m <H/R is satisfied (H is the scale-height of the 
disk at radius R). This means that the dragging of the external 
fields is always unimportant for thin disks, as V m ~ 1 is sug- 
gested, however, the advection of magnetic fields may be effi- 
cient for a geometrically th ick accretion disk with H ^R (e.g., 

Ferreira & Petrucci 2011b. A n alte rnative model was sug- 

gested bv lSpruit & U zdenskv (2005]) for advection of the ex- 
ternal field in the disk, in which turbulent diff usion is reduced 
by bu n dles of the large-scale magnetic field (Ste hle & Spruit! 
120011) . lLovelace eFatl (120091) suggested that the field can be 
efficiently advected inward based on the assumption of the 
surface layer of the accretion disk to be nonturbulent. The 
general relativistic magnetohydrodynamic (GRMHD) simu- 
lation of an accretion torus embedded in a large-scale mag- 
netic field showed that the mass is accreted mainly within the 
accretion disk, and the magnetic f ield flux is carried by the 
motions in the low-density corona dBeckwifh et al.ll2009h . 

Low mass accretion rate rh may lead to the accretio n 
flows to be advection-dominated (INaravan & Yill 1994 1 1995b . 
Advection dominated accretion flows (ADAFs) are sug- 
gested to be pres ent in low-lumi nosity active galactic nu- 
clei (AGNs) (see iNaravanl I2002T for' a review and refer- 
ences therein), FR I radio gal axies (e.g.. lGhisellini & C elotti 
1 200 U lCao & Rawlings 2004), or BL Lac objects (e.g. iCaol 
2003). The ADAF model can successfully explain most 
observational features of low-luminosity AGNs and black 
hole X-ray binaries in the l ow/hard state (e.g.. lLasota et alJ 
1996; Gammie et al.1 [1991 lOuataert et all 11999; Xu & Caol 
2009). ADAFs are hot and geometrically thick, which have 
relatively higher radial vel ocities than thin accretion disks 
(INaravan & Yill 1994 119951) . This implies that the advection 
of the external fields in ADAFs may be more efficiently than 
that in thin disks. 

In this work, we explore the advection/diffusion of the 
large-scale magnetic fields threading an ADAF, in which the 
compression of the accretion flow in the vertical direction by 
the magnetic field is properly considered. 

2. MODEL 
2.1. The structure of the ADAF 

The dynamics of a steady ADAF with large-scale magnetic 
fields are described by a set of differential equations, namely, 
the continuity equation, the radial and azimuthal momentum 
equations, and the energy equation. In this work, we use 
cylindrical coordinates (R, <fi,z). 

The continuity equation is 



dR 



(pRHv R ) = 0, 



(1) 



where mass loss rate in the winds from the ADAF is ne- 
glected. 

The radial momentum equation is 



dv R , n2 ^2, Id, , B S R B Z B H dB- 
dR K pdR HsJ 2tt£ 2nY,dR' 



(2) 



where B R is the radial component of the large-scale magnetic 
fields at the disk surface, and the pseudo-Newtonian potential 
is adopted to simulate the gravity of a non-rotating black hole 
dPaczvnskv & Wiitalll980i) . The term B r B z /2ttY, is the radial 
component of the force caused by curvature of the field lines, 

which is derived from f^ H (B z /4-iTY,)dB R /dzdzby approximat- 
ing dB R /dz — B R /H. The Keplerian angular velocity is given 



by 



Ql(R): 



GM 



(R-R^R' 



(3) 



where R g = IGMjc 1 . 

The radial magnetic fields will be sheared into azimuthal 
fields by the differential rotation of the flow, which leads to 
magnetorot ational instability (MRI) and the MRI-driven tur- 
bulence (Balbus & Hawlevl [19911, 11998). For simplicity, we 
assume that the conventional a-viscosity can describe the an- 
gular momentum transportation in the accretion flow caused 
by MRI-driven turbulence. The angular momentum equation 
is 

dtl _ V R (j-jin) , 4 s 

dR vR 1 ' K ] 

where j is the specific angular momentum of the flow, j[ n is 
the specific angular momentum of the gas swallowed by the 
black hole, and the a-viscosity, 



v = ac s H, 



(5) 



is adopted. The angular momentum equation (0]l reduces to 
an algebraic equation, 



./ = 



v R Rjh 



ac s H + v R R ' 



by assuming dfl/dR ~ —Cl/R. 
The energy equation is 

, ds 



Zv R T — =Q -QT 
R dR * * 



(6) 



(7) 



where T is the temperature and s is the entropy of the gas in 
the ADAF. Equation (0 can be re-written as 



2 dc s 1 dp fv R 



('y-l)csdR pdR 



:lR 2 H 



(j-jrn) 



(8) 



where 7 is the ratio of specific heats, and the parameter / 
describes the fraction of the dissipated energy advected in the 
flow. Substituting the continuity equation (Q~|) into Equation 
([8]), we have 



dc s 



(7-l)c s dR 



1 

R 



1 dH 
~H~dR 



1 dv R fv R 



v« dR ac 3 s R 2 H 



U-jmf- (9) 



The vertical structure of the accretion flow is significantly al- 
tered in the presence of large-scale magnetic fields. In this 
work, we calculate the ve rtical structure of the ADAF follow- 
ing the approach given in lCao & Spruij (120021) . Assuming the 
accretion flow to be isothermal in the vertical direction, we 
can calculate the vertical structure of the ADAF with 



; dp(z) 

: dz 



-piWlz- 



B R dB R B r dB z 
4^~dT + 4^~dR 



(10) 



provided the shape of the magnetic field lines in the accretion 
flow is known. An additional term containing dBj dR, which 
is the vertical component of the force caused by curvature of 
the field lines, is included, because H ~ R holds at least in 
the outer region of ADAFs. For geometrically thin accretion 
disks (H <C R), dB R /dz 3> dB 7 /dR, and the term containing 
8B z /dR is therefore neglected (ICao & Spruitll2002l) 

In principle, the field line shape is computable by solv- 
ing the radial and vertical momentum equatio ns with suitable 
boundary conditions (see lCao & Spruitll2002l for the detailed 
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discussion). For the isothermal case, an approximate analyt- 
ical expression is proposed for the shape of the field lines in 
the flow: 



x j 2 ^B R B Z | £ bz B R B : | £bz£brB R B z 



Hk 



R 



R 



R-Ri- 



%-tf+^/r 2 ) 1/2 -A(l-»£) 1/2 , {U ) 



/3 = 1- 



B R B Z 



where H is the scale height of the disk, r/[ = tanh(l), and 
the inclinat ion of the field line kq = B Z /B R at the disk sur- 
face z = H (ICao & S pruit 2002). This expression can repro- 
duce the basic features of th e Kippenhahn-Schliiter model 
(Kippenhahn & Schliiter 1957) well either for weak or strong 
field cases. The vertical distribution of the accretion disk with 
magnetic fields is not exactly the Gaussian distribution of an 
isothermal disk. In this case, we define the disk scale height 
H as p(H) = p(0)exp(-l /2), the disk scale height H can then 
be evaluated numerically with the given m agnetic field line 
shape (fTTT l. As done bv lCao & Spruit! d2002l) . we use a fitting 
formula to calculate the scale height of the disk in the rest of 
this work, 



8tt(1 -e- l l 2 )(2H+f x R) P Hn 2 K k 1 
fiR 



/4 = 1 + 



/ 3 (2ff + /!*)' 



and £br is defined as 



9Bl , , Bl 



3 

(17) 

(18) 
(19) 

(20) 



Substitute Equation ( TToT l into the energy equation (O, we fi- 
nally obtain 



H _ 1 / 4e 2 . 
R 2\R 2 n 2 K +fl 



1/2 



1 dc, _ 1 
c s dR f 4 f 5 



1- 



2c 



2f x R 



/ 3 fffi|(2ff+/ 1 tf) f 3 (2H+fiR) 



1 dfl K 

n^~dY 



1 



fiR 



-Ju 



(12) 



1 dv R 

f 5 v g dR f 5 ac\R 2 H 



(j-Jm) 



where 



Cfeg ~ fiH - f 3 f 4 H)(2H+f l R) - f\RH 
f 3 f 4 f 5 RH(2H+f l R) 

2c 2 



1 / Bl £ bz B 2 \ + ££s 

/[ " 2(l-e-i/2 )K0 l4^ff^ + 4^0|J ' (B) f 3 f 4 f 5 RHn&2H + m 



and is defined as 



where 



2c 



5 



~dR—t UR) R' 



(14) 



7-1 f 3 f 4 HQl(2H+f 1 Ry 



(21) 



(22) 



which can be calculated with the balance between the advec- 
tion and diffusion of the fields in the accretion flow (see Sec- 
tion 2.2). We find that this formula can reproduce the numer- 
ical results quite well. We note that equation (TTZb reduces to 
H = c s /57k in the absence of magnetic field. 
We define dimensionless quantities by 



In the same way, we can re-write the radial momentum equa- 
tion as 



v2 R- v l,sdv R 
v R dR 

cl(h-fs) 



2 , j 2 , c 2 s B S R B Z £ hz B 2 H 



= -Rni+^- 7 + — + 



R 3 R 2ttY> 2ttER 



R - H 
R' R 



v R 

v R = —, 

c 



~ - Si 

c 



Ms 



2c 2 . 



2fiR 



i dn K 

n^~dR~ 



[(f 2 R - f 3 H)(2H +f 1 R)- f,RH]c 2 s 



2c 



R g c R 3 g c 



f x URHttl(2H+f x R) 



- B 7 
B z =—, 
Bo 



Br 

B R = — , 

B ' 



A) 



S(/?out) 

B 2 

■■piRout)/^, 



(15) 



where Bo is the strength of the putative external imposed ho- 
mogeneous vertical magnetic fields, and R out is the outer ra- 
dius of the accretion disk. 
Differentiate Equation ( fT2l . we obtain 



2c s 



dc s 



Rfi 



dv R 



1 dH _ 

H~dR~ fTf 4 HVL\\{2H + f x R) dR f 3 f 4 v R {2H+ f x R) dR 



+ 

where 
v R s = { I 
and 



f 3 f 4 RH(2H+f l R) 
, 2/ 6 c^ 

f 3 f4fsRHS%pH+f 1 K) ' f 5 ac s R 2 H 
f 6 c 2 s [(f 2 R-f3H-f 3 f 4 H)(2H+f 1 R)-f l RH] 



fefv R 2 

+ ~ ^TTZ(J-Jm) 



f 3 f 4 f 5 RH(2H+f l R) 



fiR 



1 

"74 



2c? 



2fiR 



f 3 S%H(2H+fiK) f 3 (2H+f[R) 



1 dfl K 
n^~dR~ 



k 

f 3 f 4 (2H+f x R) f 5 

2c? 



1- 



fiR 



2c 



(f 2 R-f 3 H)(2H+f 1 R)-f 1 RH 



f 3 f 4 Htil(2H+f x R) 
It is found that Equation ( l24l > reduces to 



f 3 M2H+fiR) 



-2. 



f 3 f 4 RHnl(2H+f x R) 
where 

fi = 



f 3 f 4 RH(2H+f l R) 



(16) 



V R j 



27 
7+ 1 



1/2 



Cs- 



(23) 

1/2 

Cs, 

(24) 
(25) 

(26) 



H 



8vr(l -e-y 2 ){2H+f x R)pRn\ 



in the absence of magnetic field. We re-write the equations in 
dimensionless form as follows. 
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The radial momentum equation is 

1 dv R if . ft2 , j 2 , c~ 2 _ 2rc 2 sm v R B s R B, 



, 2rHd 2 s oat v R ^ z B 2 [(/ 2 -/ 3 H)(2ff+/i)-M]c? 



r out^out%outA) 



(/6-/5)ff 
A/5 



2c: 



,t2 



/ 3 /4tf(2ff + /i) 
2/l 



let 2h~c\ 



1 dfi K 

S7 K ^ r 



fsfA^HQKlH+fO fsMs^HfliilH+fi) 

j6C 2 s [(f2-f3H-f 3 f4H)(2H+f l )-f l H] 
f 3 f 4 f 5 rH(2H+f l ) 

fefVR ~ r 2 \ 
fsac^H J 



(27) 



where 



and 



/6 

hf A {2H+f x ) fs 



n K = 



1 



1 

/ 3 / 4 (2ff+/i). 

1 rff2 K 3r-l 



L/2 



V^rV2( r -i)' n K rfr 2r(r-l)' 
The angular momentum equation (|6]l is re-written as 

The energy equation in dimensionless form is 

2c1 2/i 



Cs, 

(28) 
(29) 



(30) 



c s flfr / 4 / 5 



1 <ffl K 
f2 K 



1- 



/1 



73/4(2//+/!) 



2c 



,~2 



/ 5 vr </r / 3 / 4 / 5 r 3 //02(2//+/ 1 ) 



, (f2-f3H-hf 4 H)(2H+f l )-f l H fv R f7 - 2 

+ = = + ^r^(j-7in) ■ 



hUfcrH(2H+fd fcaclr*H 
The functions / — /g in the dimensionless form are given by 
1 cl oul v R B%B 2 M 



h-- 



1 Hlc l,ovA V R 



(31) 
by 

(32) 



1 -e" 1/2 ffoutV*,out^ut^K-So(2H+/l) 

/ 2^ br BlB 



/3 = 1" 



1 



+ £bzB R Bz + {.bz^biB^B; 

cl,oulVR B R B z 



1 -e" 1/2 // ou tv«,outdt^l«o/3o(2//+/i) ' 
/1 



/4=1 + 



(33) 
(34) 
(35) 



/s 



2c: 



,~2 



and 



/6 



7-1 / 3 / 4 / a fffi|(2ff+/i)' 
2ci 



-2. 



/ 3 / 4 r 2 //ft 2 (2//+/ 1 ) 
The dimensionless scale height of the disk is 



H 



1 / Act 



r. 1 - 



1/2 



2 V r 2 f£ 



1 



(36) 



(37) 



(38) 



2.2. Magnetic field configuration of the ADAF 

The advection/diffusion of large-scale (poloidal) magnetic 
fields is described by 

H 

d d 4nn R [ 

-H 

(39) 

where J<p(R,Zh) is the current density at z = Zh above/below the 
mid-plane of the disk, ip(R,z) is the azimuthal component of 
the magnetic potential, and the vertical velocity v z of the flow 
is neglected. Assuming the azimuthal current distribution in 
the vertical direction to be homogeneous, we have 



Jd>(R,Zh) = 



2H ' 



(40) 



where P^(R) is the surface current density at R in the disk. As 
the work bv lLubow et al.l (119941) . the magnetic field potential 
ip(R,z) = ipd(R,z) + ip oc (R,z), where tj; d (R,z) is contributed by 
the currents in the accretion flow, and ipcc(R) = BqR/2 is the 
external imposed homogeneous vertical field, which can be 
regarded as being contributed by the currents at infinity (see 
iLubow etaUl 1 9941 for the details). The potential tp& is related 
to PAR) with 



2tt 



1 



^ d (/?,z) = - / R'dR' j cos W 

«in 



J*(R',Zh) 



[R' 2 + R 2 + (z-z h ) 2 - 2RR' cos </>'] 1 I 2 



— j P^R')R'6R' j cos W 



dzh- 



[R' 2 + R 2 + (z - z h ) 2 - 2RR' cos </>'] l l 2 
—ti 

Differentiate Eq. (|4TT >. we have 

Rout 2-K 

^[^d(^,z)] = ^ / PAR')R'dR' J cos0'd<// 

Kin 



dzh- (41) 



R'"- + (z-Zh) 2 -RR' cos 0' 
[R' 2 + R 2 + (z-z h ) 2 - 2RR> cos 0'] 3 /2 



dz h . (42) 
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For steady case, i.e., d/dt = 0, equation (f39l > becomes 



d r , , 2ir ac s R r 

■j-[RMRM 

OK c V R 



r m 4(R) = B R, (43) 



where the magnetic Prandtl number is define as V m = r\jv. 
This equation can reduce to a set of linear algebraic equations, 
i.e., 

-j2Pu4(R J )ARj--^^V m 4(R i ) = B oRi , (44) 



7=1 



v R (Rd 



by using Equation d42l) . where J^(Rj) is the surface current 
density of the ring at radius Rj in the accretion flow, ARj is 
the width of the ring, and 



Pir 



2H jC 



COS( 



R 2 + z h 2 



- RjRj cos < 



[Ri 2 + Rj 2 + zl~ 2RiRj cos </>'] 3 /2 



dzh- 



(45) 

The surface current density J^{Rd can be calculated by solv- 
ing a set of linear algebraic equations provided the structure of 
the ADAF i s given with the sp ecified magnetic Prandtl num- 
ber V m (see lLubow et al-lfl 9941 for the details), and therefore 
the configuration of the large-scale magnetic fields is available 
with the derived potential ip: 



and 



B R (R,z) = ~ip(R,z), 
oz 



B z (R,z)=~-^m(R,z)l 
2.3. Boundary conditions 



(46) 



(47) 



In order to carry out the calculation of the structure of the 
ADAF, we have to specify the boundary conditions at the 
outer radius R out . There are six model parameters: viscos- 
ity parameter a, the outer radius of the accretion disk r out , the 
de gree of advection _/", specific angular momentum j at r out , 
the temperature of the gas 8 out = c s /Rfl^ = V2(r- l)c s /r 1//2 
at r out , and /?o = %Kp(R out )/B\, which relates the external im- 
posed vertical magnetic field strength to the gas pressure of 
the disk at r out . In all the calculations, we fix the ADAF pa- 
rameters: a = 0.2, / = 0.95, j(r out ) = jK(r ut), and 9 out = 0.25, 
which are typical for ADAFs. We find that the final results 
are insensitive to the values of these parameters. We integrate 
Equations d27l ) and ( f3Tb inwards from R = R oul numerically. 
The parameter j m , the specific angular momentum of the gas 
swallowed by the black hole, is tuned carefully till the derived 
global solution passes the sonic point smoothly. 

The dynamical structure of the ADAF is coupled with the 
magnetic field, i.e., the structure of the ADAF is affected by 
the large-scale field threading the flow, and vice versa. In 
principle, the derived solution has to pass through several crit- 
ical points (e.g., the slow/fast magneto-sonic points and the 
Alfven point). In this work, we first calculate the structure 
of the ADAF without magnetic field, and the magnetic field 
configuration/strength is then calculated based on the derived 
ADAF structure. The structure of the ADAF is re-calculated 
with this derived field configuration/strength. The calcula- 
tions of the dynamical structure of the ADAF is decoupled 
with those of the field configuration, and the global solution 
of the ADAF is therefore only required to pass the sonic point. 



The calculations are iterated till the solutions converge. We 
find that the final solution is available usually after three or 
four iterations. This is because the dynamics of the ADAF is 
altered little by the magnetic field except in the inner edge of 
the flow. 



3. RESULTS 

The configuration of the large-scale magnetic fields are 
plotted in Figures Q] and [2] The structure of ADAFs is al- 
tered due to the presence of such large-scale magnetic fields, 
which is plotted in Figures |3]|6] It is found that the radial 
velocity of the accretion flow changes little except in the re- 
gion of the accretion flow close to the black hole horizon (see 
Figures [3] and H}. The accretion flow is decelerated near the 
black hole by the magnetic field when f3o is relatively small 
(see Figures j3][6]l, which implies that the accretion flow can 
be trapped by the magnetic field provided the external im- 
posed field is strong enough. We find that the accretion flow 
may be disrupted by the magnetic field when the parameter 
/3o is lower than a critical value /3o.cnt, provided all other disk 
parameters are fixed. The critical values of (3q as functions of 
outer radius of the ADAF are plotted in Figure [7] The verti- 
cal structure of the accretion flow is altered in the presence of 
magnetic fields, i.e., the scale height decreases significantly in 
the inner region of the accretion flow (see the middle panels 
in Figures |5Jj6]l. In Figures I8TTPT1 we compare the relative im- 
portance of the gas pressure and the magnetic pressure in the 
accretion flow with different magnetic field parameters (i.e., 
/3o and V m ). The ratios of the magnetic field strength in the 
accretion flows to the external imposed field strength as func- 
tions of radius are plotted in Figure [121 
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FIG. 1 . — The large-scale poloidal magnetic field configuration of an ADAF 
(the magnetic field strength 0g = 200 is adopted). The dotted line is the 
scale height of the ADAF. The outer radius of the disk is assumed to be 
Rout = lOORg. The magnetic field configuration plotted in the upper panel 
is calculated for the magnetic Prandtl number V m = 1 , while the lower panel 
is for Vm = 1.5. Every magnetic field line corresponds to a certain value of the 
stream function Rip(R,z) with the lowest [Rip(R,z)] m m =0.04R ollt i/ , oo(Sout,0) 
increasing by ARi(>(R,z) = 0.04K out i/>ooCKout,o) P er une . where ^oo(Rout) = 





FIG. 3. — The structure of the accretion flows. The magnetic Prandtl num- 
ber Vm = 1 is adopted in all the calculations. The color lines represent the 
results with different magnetic field strengths, i.e., /3q = 20 (red), 40 (green), 
200 (blue), and 400 (yellow). The black lines represent the structure of the 
accretion flow without magnetic fields. The upper panel: the radial velocities 
(solid lines), and the sound speeds (dotted lines) as functions of radius. The 
middle panel: the scale heights of accretion flows as functions of radius. The 
lower panel: the specific angular momenta as functions of radius. The black 
dashed line represents the Keplerian specific angular momentum. 



The large scale mag: 




FIG. 4. — The same as Figure [5] but the outer radius of the accretion flow 
Rout = 1000/fg is adopted. The color lines represent the results with different 
magnetic field strengths, i.e., j3o = 2 (red), 4 (green), 20 (blue), and 40 (yel- 
low). 
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FIG. 5. — Similar to Figure [5] the results with different values of the mag- 
netic Prandtl number are compared (/3q = 200 is adopted). The red lines rep- 
resent the results calculated with V m = 1 , while the green lines are for those 
with Vm = 1.5. 



4. DISCUSSION 

The configurations of the magnetic fields show that the 
advection of magnetic fields is ve ry efficient if V m ~ 1 is 
adopted, unlike the thin disk cases dLubow et alj [l994). This 
is because the radial velocity of ADAFs is much higher than 
that of the thin disk. Such configurations may help launch- 
ing outflows/jets from ADAFs. For simplicity, we have not 
considered magnetically driven outflows from ADAFs in this 
work. The radial velocity of the ADAF will increase if the 
angular momentum carried away by the outflows is properly 
taken into account, which will enhance the advection of the 
magnetic fields in the ADAF. We have not considered the az- 
imuthal component of magnetic field B<p in the radial momen- 
tum equation (O, which may affect the dynamics of the accre- 
tion flow. Our results show that the dynamics of the accretion 
flow is dominantly determined by the gas pressure term ex- 
cept in the region very close to the black hole horizon. The 
previous MHD simulations showed that the magnetic stress 
never exceeds, usually much lower than, t he gas pressure in 
the mid-plane of the accretion flow (e.g., [Hawley & Krolikl 
1200 U |i gumenshchev et al. 2003), which implies that the main 
conclusions of this work will not be altered even if the az- 
imuthal field component is included in the radial momentum 
equation. 
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FIG. 6. — The same as Figure [5] but the outer radius of the accretion flow 
Rout = lOOORg, and the field strength parameter fig = 20 are adopted. 




FIG. 7. — The critical values of @q as functions of outer radii of ADAFs 
with Vm = 1 (red) and 1.5 (green), respectively. 




FIG. 8. — The magnetic pressure/gas pressure of the accretion flow, and the 
inclination of field lines at the disk surface as functions of radius. The mag- 
netic Prandtl number Vm = 1 is adopted in all the calculations. The different 
color lines represent the results with different magnetic field strength, i.e., 
00 = 20 (red), 40 (green), 200 (blue), and 400 (yellow). The upper panel: 
the magnetic field pressure/gas pressure of the accretion flow as functions 
of radius. The black line represents the gas pressure of the ADAF without 
magnetic fields as a function of radius. The middle panel: the inclination of 
the magnetic field lines at the disk scale height. The lower panel: the ratio 
of magnetic pressure to the gas pressure in the accretion flow (solid lines), 
and the ratio of magnetic pressure to the gas pressure in the accretion flow 
without magnetic fields (dashed lines). 

The ratio of the magnetic pressure to gas pressure can reach 
~ 0.5 in the inner edge of the ADAF (see Figures l8lfTTTi. The 
radial velocity and the temperature distribution have been al- 
tered little in the presence of large scale magnetic fields, ex- 
cept in the inner region of the ADAF (< 3R g ). However, the 
ADAF is vertically pressured by the magnetic fields, which 
leads to the gas pressure significantly increased in the inner re- 
gion of the ADAF. It is found that the magnetic pressure of the 
ADAF in the region close to the black hole horizon can be two 
orders of magnitude higher than the gas pressure of the ADAF 
without magnetic fields (see the dashed lines in the lower pan- 
els in Figures [8] and [9). In most previous works, the mag- 
netic field strength is estimated with the gas/radiation pressure 
in the accretion disk on the so-called equipartition assump- 
tion, and the confinement of the disk by the magnetic fields 
in th e vertical direction of the flow has not been considered 
(e.g.,|M oderski & Sikorall 19961: iGhosh & Abramowiczll 19971 
Livio et alJ1999t | Armitage & N ataraian 1999; Nemme n et al.l 
2007tlWu & Caoll2008tlMcNamara et alJl2009b . It means that 
the maximal jet power estimated based on the normal ac- 
cretion disk models without considering magnetically con- 
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finement in vertical direction is around two orders of magni- 
tude underestimated. The assumption that the strength of the 
magnetic fields near the black hole horizon should not dif- 
fer significantly in strength with th e accretion disk m ay be 
incorrect, at least for ADAFs (e.g. iLivio et ail ll 999). Our 
present calculations are based on the pseudo-Newtonian po- 
tential, which can only simulate the gravitational potential of 
non-rotating black holes. However, the gas in the accretion 
flow plunges rapidly to the black hole in the region of the 
flow with R < R ms (R ms is the radius of the marginal stable 
circular orbit of the hole) for a rotating black hole, which is 
similar to the cases calculated in this work for non-rotating 
black holes. We believe the conclusion that the magnetic 
fields are significantly strengthened near the black hole hori- 
zon due to very large radial velocity of the gas in the flow 
should still hold even for rapidly spinning black holes. This 
implies that the efficiency of the Blandford-Znajek mecha- 
nism is sig nificantly underestimated in most of the previous 
work s (e.g..lModerski & S ikora 1996; Ghosh & Abra mowicz l 



cz 
ao 



1997; Livi oetalJ 11 999; Ne mmen et all I2007t IWu & Cat. 
20081; iMcNamara et al.ll2009l) . The detailed calculations for 
rotating black holes is beyond the scope of this work, which 
will be reported in our future work. 
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FIG. 10. — Similar to Figure[8] the results with different values of the mag- 
netic Prandtl number are compared (/?o = 200 is adopted). The red lines rep- 
resent the results calculated with Vm = 1 , while the green lines are for those 
with Vm= 1.5. 



The scale height of the ADAF is sensitive to the external im- 
posed homogeneous vertical magnetic field strength. We find 
that the scale height of the ADAF decreases with increasing 
field strength, i.e., low-/3o case (see the middle panels in Fig- 
ures [3] and |4j. The scale height also depends on the value of 
the magnetic Prandtl number V m . The scale height increases 
with the value of V m (see Figures [5] and |6j, because magnetic 
fields being less stronger for a higher V m (see Figures [To] and 
ITTb . The relative thickness of the disk H/R becomes very 
small (i.e., geometrically thin) in the inner edge of the ADAF, 
due to strong magnetic pressure. It implies that, the optical 
depth in the vertical direction will increase significantly, and 
the temperature difference between ions and electrons will de- 
crease in the inner edge of the ADAF, though no detailed ra- 
diative processes are included in our present calculations. Its 
implications on the observational features of the objects ac- 
creting at low rates can be explored, if the detailed radiative 
processes are considered and the black hole mass and the mass 
accretion rate need to be specified, which is beyond the scope 
of this work. 



FIG. 9. — The same as Figure [I] but the outer radius of the accretion flow 
Rout = lOOOSg is adopted. The different color lines represent the results with 
different magnetic field strength, i.e., /3q = 2 (red), 4 (green), 20 (blue), and 
40 (yellow). 
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FIG. 1 1 . — The same as Figure [Tol but the outer radius of the accretion flow 
Rout = lOOORg, and the magnetic field strength parameter /?□ = 20 are adopted. 




Rl R 



FIG. 12. — The ratios of the magnetic field strength in the accretion flows to 
the external imposed field strength as functions of radius. The red lines repre- 
sent the results with /3q = 20 for the accretion flows with R out = iOORa (/3o = 2 
for the accretion flows with R out = 1000K g ), while the blue lines represent the 
results with 0o = 200 for the accretion flows with R out = 100R g 0o = 20 for 
the accretion flows with R out = 1000S g ) The solid lines represent the results 
for V m = 1, while the dashed lines are for V m = 1.5. 

It is interesting to find that the accretion flow is decelerated 
near the black hole by the magnetic field when /3o is relatively 
small (see Figures [3] and |4). There is a critical value /3o,crit, 
below which the accretion flow is decelerated near the black 



hole. The value of /3o,cnt decreases with increasing outer ra- 
dius R out , and it decreases with increasing the value of the 
magnetic Prandtl number V m (see Figure|7]i. This implies that 
the accretion flow may be disrupted by the magnetic field in 
the inner region of the accretion flow, when the external im- 
posed field is strong enough or the gas pressure of the flow is 
low at the outer radius, or both, which justifies the main as- 
sumption in the qualitative analysis o n the magnetic arrested 
accretion disks (Naravan et al. 2003). In this case, the gas 
may accrete as magnetically confined blobs diffusing through 
field lines in the region v ery close to the black hole horizon 
(see iNaravan et al.H2.Q03l for the detailed discussion, and the 
refer ences therein), which is similar to those in compact stars 
(e.g jElsner & Lambll984UKaisig et alll992h . For an ADAF 
surrounding a black hole with a given external imposed ho- 
mogeneous field, the value of /3o,crft corresponds to a certain 
accretion rate at R out , below which the accretion flow will be 
trapped by the magnetic field, for given temperature of the 
gas in the flow at R out . Our calculations show that the gas 
can alternatively be trapped by the magnetic field at outer 
radius of the accretion flow provided the magnetic pressure 
significantly exceeds over the gas pressure, i.e., (3q <C 1 (see 
the green line in Figure [7] at large radii). We compare the 
amplification of the magnetic field advected in the accretion 
flows with different values of the disk parameters, and find 
that the amplification of the field is dominantly determined 
by the value of V m , while it is insensitive to the adopted value 
of (3q. It implies that the strength of the field near the black 
hole horizon is mainly determined by the strength of the ex- 
ternal imposed field and the outer radius of the ADAF, which 
is quali tatively consistent with the results for thin accretion 
disks in lLubow et alj (119941) . though the advection of the field 
is rather inefficient in the thin disk case. This means that the 
strength of the magnetic field near the black hole can be esti- 
mated if the strength of the ordered field threading the ambi- 
ent gas and the outer radius R out of the ADAF are known. It is 
found that stronger magnetic field will be in the inner edge of 
the ADAF, if the gas starts to accrete from a larger radius (see 
Figure [T21. This should be very useful for estimate of the field 
strength near the black hole horizon in some low-luminosity 
active galactic nuclei, of which the properties of the ambient 
gas at the Bo ndi radius have well been measured w ith obser- 
vations (e.g. jBaganoff et alJl2003l:lAllen et al.ll2006h . 

We notice that a series of work s of numerical simulations 
on th e magnetic fields in A DAF s dlgumenshchev & Naravanl 
tlOOH flgumens hchev et al.ll2003Ulgumenshchevll2006r i2008). 
They found that the initial uniform z-direction magnetic field 
can be efficiently dragged inwards by the ADAF, and their de- 
rived magnetic field config uration is similar to those obtaine d 
in this work (see Figure 3 in Igumenshc hev & Naravanl2002|) . 
The numerical simulations show that the magnetic field in the 
inner edge of t he accretion can be ve ry strong to disrupt the 
accretion flow (Igumenshchev 2008), which agrees qualita- 
tively with our calculations. Compared with these MHD sim- 
ulations on ADAFs, our calculations explicitly indicate that 
the vertical scale-height of the ADAF near the black hole 
is significantly reduced, and hence the field strength can be 
much higher than that estimated based on the conventional 
ADAF model. 

In this work, the viscosity a = 0.2 is fixed. The radial ve- 
locity of the ADAF in the outer region is roughly proportional 
to the kinetic viscosity v, and the magnetic Prandtl number 
V m = rj/v, which implies the advection/diffusion process in 
the accretion flow is almost independent of the adopted value 
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of a. The dynamical structure of ADAFs is insensitive to the 
value of / provided the accretion flow is advection dominated, 
i.e., 0.5 < / < 1. We also calculate the problem with different 
values of the disk parameters (e.g., /, out , and j out ), and find 
that the main conclusions have not been altered. 

The ADAF may connect to a thin accretion disc at a 
certain transition radius R tI . This is required by model- 
ing on a variety of observations of AGNs/X-ray binaries 



IU997tlOua taert et a fll 19991: iLu & Wang|l2000t 
ICaol I2003h lYuan & Narayanll2004t IXu & Cad 120091 and is 
also predicted by some theoretical model calculations (e.g., 
Abram owicz et all 1 9951 iLiu et al.ll999URozaNska & Czernvl 
2000tlSpruit & DeufeJ2002ElLu et al. 2004). Our calculations 



are also valid for the case that the inner ADAF connects to the 



outer cold thin disk at a certain radius. In this case, the advec- 
tion of the external fields is quite inefficiently in the outer thin 
disk due to its low radial velocity, and the field lines thread 
the disk almost vertically (see, e.g. JLubow et al.| [l994), while 
these field lines can be efficiently dragged inward by the radial 
motion of the inner ADAF. 
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